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^ ! 1 INTRODUCTION 

b> ' 

, J_| ■ Weak gravitational lensing by the cosmic-large scale structure pro- 
' vides a measurement of the cosmic tidal field and provides sensitiv- 
, ity on cosmological parameters due its dependence on the geometry 
&3 • of the cosmological model and the growth rate of matter pertur- 
bations (for reviews, see Bartelmann & Schneider 2001; Hoekstra 
& Jain 2008; Bartelmann 2010). Weak lensing data taken by fu- 
ture surveys such as EUCLID, DES and LSST has the potential of 
putting tight constraints on cosmological parameters and to distin- 
guish between dark energy models and those with modified gravity 
(for a comprehensive summary, see Weinberg et al. 2012). 

In analysing weak lensing data it is usually assumed that weak 
lensing is the only effect causing correlations in the shapes of galax- 
ies, which is sadly not the case as correlations in the shapes of 
neighbouring galaxies are naturally explained by their correlated 
angular momenta due to similarities in their formation processes, 
most importantly the tidal shearing experienced by their host haloes 
(for a review on angular momentum models and intrinsic align- 
ments, see Schafer 2009). 

Angular momentum is introduced into haloes in their forma- 
tion as a consequence of tidal shearing which has been investigated 
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ABSTRACT 

Subject of this paper is the weak lensing effect on galaxies that show intrinsically correlated 
ellipticities. In our model, we investigate the distortion of the ellipticity field if the galaxies 
experience an apparent shift in their position by weak lensing deflection and compare this 
effect to the shearing effect induced by tidal fields. Starting with a derivation of intrinsic ellip- 
ticity spectra by employing a tidal torquing model generating galactic angular momenta, we 
model the galaxy ellipticity by assuming that the galactic disk forms perpendicularly to the 
host halo angular momentum direction and derive intrinsic ellipticity E-mode and B-mode 
spectra from the angular momentum statistics. The lensing effect on the ellipticity field is 
modeled by employing the methodology developed in the framework of lensing of the cosmic 
microwave background polarisation. For EUCLID, ellipticity correlations are altered by lens- 
ing deflection on multipoles i i 1000 by ~ 5% for the ellipticity Zs-modes and by ~ 30% for 
the B-modes, while a shallower survey would exhibit larger changes on larger angular scales. 
In addition to the convolving effect of lensing on the ellipticity spectra we investigate the 
E/B-mode conversion, and discuss the possibility of measuring correlations between different 
multipoles which is evoked by the homogeneity breaking effect of the lensing displacement. 
Our conclusion is that although shape correlations generated by weak gravitational shear is 
dominant, the shifting effect due to lensing is shaping the ellipticity spectra on small angular 
scales and causes a number of interesting phenomena, which might be observable by future 
surveys. 

Key words: cosmology: large-scale structure, gravitational lensing, methods: analytical 



using perturbation theory, including the deformation of forming 
haloes, as well as by numerical studies (Catelan 1995; Catelan & 
Theuns 1996b,a; Catelan et al. 2001; Catelan & Porciani 2001; Lee 
et al. 2007; Codis et al. 2012). There, the role of tidal torquing in the 
angular momentum build-up during galaxy formation is supported, 
but indicate that tidal torquing models might be predicting too high 
values for the amount of angular momentum (Bullock et al. 2001; 
Catelan & Porciani 2001 ; Porciani et al. 2002a,b; Hahn et al. 2007) 
while perturbation theory seems to yield fairly reliable results for 
the angular momentum direction and explaining alignment effects 
of the halo with the ambient large-scale tidal field. 

Angular momentum alignments give rise to a correlation in 
intrinsic ellipticities between neighbouring galaxies if the angular 
momentum direction of the galaxy corresponds to the one of the 
host halo. The scale of this correlation is predicted to be in the 
range of ~ 1 Mpc/h (Crittenden et al. 2001; Natarajan et al. 2001; 
Schafer & Merkel 2012). Hence, intrinsic alignments influence the 
angular momentum correlation on small scales and can be of sig- 
nificance for high precision observations such as the future surveys 
mentioned above. Estimates of intrinsic alignment spectra arising 
from angular momentum based-models can be found in the works 
of Croft & Metzler (2000); Crittenden et al. (2001, 2002); Mackey 
et al. (2002). Due to their large complexity, Schneider & Bridle 



©2013RAS 



2 A. Giahi-Saravani, B.M. Schafer 



(2010) devised a model which allows a much easier handling of 
intrinsic alignments. 

Direct identification of the symmetry axis of the disc with 
the host halo angular momentum would lead to overestimation of 
the ellipticity alignments and is flawed as there might be large 
misalignments between these two vectors due to baryonic physics 
(Bailin et al. 2005; Stewart et al. 2013). Thus, we want to think of 
our ellipticity spectra as upper limits. Whereas the alignment model 
for spiral galaxies a quadratic dependency of the ellipticity on the 
tidal field is assumed Crittenden et al. (2001), elliptical galaxies fol- 
low a simpler linear alignment mechanism (Hirata & Seljak 2004). 
Sadly, the processes determining the orientation of the stellar disk 
inside a dark matter halo are not amenable to analytic calculations 
and even numerical simulations struggle to reach consensus about 
properties of galactic disks in a cosmological volume. 

Intrinsic alignments are difficult to quantify (Godlowski 2012) 
but detections have been reported in a number of data sets, for in- 
stance in the Tully-catalogue (Pen et al. 2000), the Point Source 
Catalogue Redshift survey (Lee & Pen 2002), in the Sloan Digi- 
tal Sky survey (Hirata et al. 2004; Mandelbaum et al. 2006; Lee & 
Pen 2007; Okumura & Jing 2009; Joachimi et al. 201 1), in com- 
bination with the 2-degree Field Galaxy Redshift survey (Hirata 
et al. 2007) and the WiggleZ Dark Energy Survey (Mandelbaum 
et al. 201 1), which allowed tests of alignment models (Blazek et al. 
201 1) and the determination of alignment model parameters. Man- 
delbaum et al. (201 1) in particular provide a detailed description of 
the technical difficulties in measuring intrinsic correlations at red- 
shifts relevant for gravitational lensing, and give fits for the intrinsic 
ellipticity correlation and the ellipticity-density cross-correlation. 
Some studies, however, remain sceptical about these detection re- 
ports (Andrae & Jahnke 201 1). 

As observations of galaxy ellipticities become increasingly 
precise, the study intrinsic alignments of galaxies based on angular 
momentum models will gain much interest in the future. Such mod- 
els (Croft & Metzler 2000; Crittenden et al. 2001, 2002; Mackey 
et al. 2002) have been applied to estimate the significance of intrin- 
sic alignments in the weak lensing convergence spectrum (Heav- 
ens et al. 2000; Heymans & Heavens 2003; Heymans et al. 2004) 
and bispectrum (Semboloni et al. 2008). Techniques for separating 
weak lensing data from intrinsic alignments range from discard- 
ing close galaxy pairs (King & Schneider 2002, 2003) to introduc- 
ing a weighting in order to null out their contribution (Joachimi & 
Schneider 2008; Kitching & Taylor 201 1), amplifying them relative 
to the weak lensing-induced correlations (King & Schneider 2003; 
Joachimi & Schneider 2010a,b; Joachimi & Bridle 2010) or to use 
the non-zero vortical modes (B-modes) of the intrinsic ellipticity 
field (Crittenden et al. 2002). With a model for intrinsic alignments 
at hand, it is of course possible to fit both the weak lensing elliptic- 
ity correlations and the intrinsic correlations at the same time, pos- 
sibly constraining intrinsic alignment parameters (Bernstein 2009; 
Kirk et al. 201 1) or use self-calibration techniques (Zhang 2010). 
Many of the removal techniques require very good control of red- 
shift estimates (Bridle & King 2007) and can be extended to deal 
with cross-correlations between weak lensing and intrinsic align- 
ments (King 2005). A puzzling result is that the contaminating ef- 
fect of intrinsic alignments predicted from different angular mo- 
mentum models in weak lensing data can be dramatically differ- 
ent. If intrinsic alignments are present in weak lensing data but not 
included in the model for the data when determining cosmolog- 
ical parameters, biases in the parameter estimates have to occur: 
Whereas Kirk et al. (2010, 2012) using a linear model found most 



biases in the dark energy equation of state, Capranico et al. (2012) 
arrive at the conclusion that biases arise mostly in n,„ and cr 8 . 

Motivation for this paper is to study the influence of weak 
gravitational lensing on the observed ellipticity correlation, aris- 
ing from the lensing effect of the intervening matter distribution. 
In order to assess the significance of intrinsic alignments in weak 
lensing measurements, Crittenden et al. (2001) found that the in- 
trinsic signal is between 1 to 10 per cent of the measured lens- 
ing signal for a deep reaching survey. Here, our aim is to inves- 
tigate the implications of the distortion (caused by weak gravita- 
tional shear) of galaxy ellipticities and the shift of galaxy positions 
by weak lensing on the shape of the E- and B-mode spectra of 
the lensed ellipticity correlation function. We will thus employ a 
formalism based on lensing of the cosmic microwave background 
(CMB) polarization (Seljak 1996; Hu & White 2001; Challinor & 
Lewis 2005; Lewis & Challinor 2006). It should be mentioned that 
in contrast to the lensing of the CMB where only the shifting effect 
is of importance and the polarization of the CMB photons remains 
unchanged, galaxy lensing involves also a distortion of ellipticities 
and a change of their orientation. In this regard, the analogy be- 
tween the polarization of the CMB and the orientation of galaxy 
ellipticities has to be viewed differently. We use the model of Crit- 
tenden et al. (2001) for deriving a correlation function for the in- 
trinsic ellipticities and apply the CMB formalism to it in order to 
observe characteristic features of weak lensing, such as the mixing 
of E- and B-modes. We expect a suppression of ellipticity corre- 
lations due to the shifting effect which dilutes and randomises the 
galaxy position at which ellipticities are measured. 

In his paper we collect necessary results from cosmology, 
angular momentum generation in galaxies, ellipticity correlations, 
and weak lensing in Sect. 2, before describing in detail the formal- 
ism we use for computing lensed ellipticity spectra in Sect. 3, along 
with an investigation of the effects predicted by our formalisms for 
the case of a high- and a low-redshift galaxy sample. We summarise 
our main results in Sect. 4. As reference model we work with a ba- 
sic spatially flat vvCDM model with Gaussian adiabatic initial per- 
turbations in cosmic density field. Model parameters were set to be 
n,„ = 0.25, n s = 1, erg = 0.85, Q. b = 0.04 and the Hubble-radius 
c/H = 2996.9 Gpc/h, with h = 0.72. The dark energy equation of 
state parameter w was assumed to be constant at a value of -0.9. 



2 COSMOLOGY 

2.1 Dark energy cosmologies 

The time evolution of isotropic Friedmann-universe with homoge- 
neous dark matter and dark energy is described by the Hubble func- 
tion H(a) = dlna/d/, which is given by 



— — = — +(1 -fi,„)exp 3 dlna 



(l+w(a))\, 



(1) 



with the matter density parameter Q m and the dark energy equation 
of state function w(a). Spatial flatness requires the dark energy den- 
sity to be 1 - fl,„. The comoving distance^ can be computed from 
the scale factor a, 



f ' da 

X a 2 H(aY 



(2) 



For the galaxy redshift distribution n(z)dz, we use a standard shape 



n(z) = no I — | exp 



dz with - = ^r(l , ( 3) 

«o P \PI 
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with /? = 3/2. We choose zo such that the distribution has a me- 
dian redshift of 0.9 corresponding to EUCLID (Amara & Refregier 
2007) which we contrast with a galaxy distribution of identical 
shape but with a much lower median of 0.3. We will refer to the 
two application cases as the high and low redshift galaxy sample, 
respectively. The distribution can be rewritten in terms of comoving 
distance using the relation p(z)dz = p(x)dx w i m dz/dx = H(x)/c. 

2.2 CDM power spectrum 

The statistical properties of the overdensity field 6 is in the case of 
homogeneous, isotropic and Gaussian fluctuations described by the 
spectrum Ps(k). Inflationary models suggest the ansatz 

P s (k) oc k'°T 2 (k), 



(4) 

with the transfer function T{k). This function is approximated by 

(5) 



ln(l + 2.347) i 
r(g)= 2.34, M ^> 
(see Bardeen et al. 1986), with the polynomial 
M(q) = 1 + 3.89<? + (16.1?) 2 + (5A6qf + (6.71?) 4 



(6) 



The transfer function depends most strongly on the parameters Ci„ 
and h which form the shape-parameter T, 



T = Q,,„h exp -Q, 



1 + 



(7) 



and shows slight corrections due to the baryon density Q. b 
Sugiyama (1995). With T, the wave vector k is scaled according 
to q = k/T. The linearly evolved spectrum P(k) is normalised to 
show the variance o"g on a comoving scale of R = 8 Mpc/h, 



f tfdk 



P s (k)W 2 (kR) 



(8) 



with a Fourier transformed spherical top hat filter function, W(x) - 
3 ji(x)/x. jc{x) is the spherical Bessel function of the first kind of 
order I (Abramowitz & Stegun 1972). Weak lensing spectra will 
be computed for both linear and nonlinear structures, for the latter 
we employ the extension of P(k) to nonlinear structure formation 
proposed by Smith et al. (2003). In the limit of small perturbations 
|i5| <k 1, cosmic structure formation is linear and homogeneous, 
S(x,a) = D + (a)6(x,a = 1). The time evolution of the cosmic den- 
sity field is then given by the growth function D + (a), that follows 
from the growth equation (Turner & White 1997; Wang & Stein- 
hardt 1998; Linder & Jenkins 2003), 



d 2 1 , 

— D + (a)+ - 3 + 
da l a 



^-^) ^-D + (a) = ^-a m {a)D + (a\ 
dm a / da 2a z 



(9) 



such that the spectrum P(k) oc D 2 (a), while the nonlinear extension 
to P(k) has its own parametrised time evolution based on £l„,(a), 

1 Hi 



Q.„,{a) 



Q m a 3 H 2 (a)' 
which can be derived from the adiabaticity of the cosmic fluids. 



(10) 



, dD + 
L„ = a H(a)——e a p Y Ip S <b# y , 



(11) 



with the tidal shear being defined as the Hessian of the gravitational 
potential and the inertia tensor measuring the second moments of 
the mass distribution inside the protohalo, 



<9 2 <D 



and 



dx 6 dx y 

Ips = Q,„Pc„t J d 3 q6(q)(q- q)p(q - qh, 



(12) 
(13) 



respectively, with implicit summation over repeated indices, q are 
Lagrangian coordinates moving along with the halo's centre of 
gravity q. Note that the second moments of the mass distribution 
Ips is referred to as the inertia tensor despite the fact that it has a 
different ordering of the axes. 

This relation reflects the interesting requirement of misalign- 
ments between the shear and inertia eigensystems which is nec- 
essary for angular momentum generation: Only the antisymmetric 
tensor X^ y = SaC/M^y ~ ® fishy) is relevant for the angular mo- 
mentum (Schafer & Merkel 2012), L a oc X^, because the contrac- 
tion of the symmetric contribution Xt, with the antisymmetric e a ps 
vanishes. The antisymmetric tensor X- is equal to the commutator 
[Ips, O.sy] which suggests that the tidal shear and the inertia are not 
allowed to be simultaneously diagonalisable and must not have a 
common eigensystem, otherwise angular momentum can not arise. 

In this work we use the angular momentum-based ellipticity 
correlation model proposed by Crittenden et al. (2001). There, el- 
lipticities are set into relation to tidal shear by means of a condi- 
tional probability distribution p{L\<b „p)dL. Such a distribution has 
been proposed by Lee & Pen (2001) as being Gaussian with the 
covariance 



cov(L)„« oc - 



1 /l +fl 



■ a (i> 2 ), 



up 



(14) 



which acquires a dependence on the tidal shear tensor, that has been 
normalised and made trace-free, tr(O) = and tr(0 2 ) = 1. In this 
way, the variance of the angular momentum field varies with the 
tidal shear, and the randomness of the angular momentum field 
is controlled by the misalignment parameter a, which describes 
the average orientation of the protohalo's inertia to the tidal shear 
eigensystem. a has been measured in numerical simulation to be 
the value = 0.25, and we will use this value in this work. 

This description is valid on scales where the correlations be- 
tween inertia tensors are negligible and is sufficient because we are 
only interested in the angular momentum direction, as will be ex- 
plained in the next section. Therefore, one does not need the vari- 
ance (L 2 ) of the angular momentum field as a parameter and it is 
possible to marginalise the distribution over the magnitude L, 



p(M®. a p) 



J 



L dL p(L\®, c 



(15) 



2.3 Angular momentum from tidal shearing 

Angular momentum generation in CDM-haloes is a Lagrangian 
perturbative process (Hoyle 1949; Sciama 1955; Peebles 1969; 
Doroshkevich 1970; White 1984): The variation of velocities dis- 
placing a protogalactic region acts as a torque which generates an- 
gular momentum L a : 



In this picture correlations between angular momenta can be traced 
back to correlations between tidal shears that neighbouring galax- 
ies experienced in building up their angular momenta. It should be 
noted, however, that recent investigations provide evidence that the 
nonlinear evolution of galactic angular momenta can be vorticity 
driven (Libeskind et al. 2012; Lee 2013) as an alternative to tidal 
torquing. 
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2.4 Intrinsic ellipticity correlations 

We assume that the halo angular momentum axis and the symme- 
try axis of the galactic disk are parallel. The tilting of the disk rel- 
ative to the line of sight determines the aspect ratio under which 
the galactic disk is viewed and therefore the galaxy's ellipticity. In 
this way, the angular momentum direction L = L/L with L = \L\ 
determines the ellipticity e. In this picture, ellipticity correlations 
are derived from angular momentum direction correlations and ul- 
timately from tidal shear correlations (Heavens et al. 2000; Critten- 
den et al. 2001, 2002; Mackey et al. 2002; Heymans & Heavens 
2003). Specifically, the two components of the ellipticity field can 
be combined to form the complex ellipticity with e + as the real and 
e x as the complex part, and related to the angular momentum direc- 
tion L: 



with 



1 +L? 



2a 



L X Ly 

l + £?' 



(16) 



when the coordinate system is aligned with its z-axis parallel to the 
line of sight. Under a rotation of the coordinate frame by an an- 
gle <p the complex ellipticity transforms according to the relation 
e — > exp(2i</j)e, in accordance with the spin-2 property of the el- 
lipticity field, a is the disk thickness parameter and can be used 
for describing a weaker dependence of e on L if the galactic disk 
has a finite thickness. We use the numerical value a = 0.75 which 
has been measured in the APM galaxy sample by Crittenden et al. 
(2001). 

Linking the ellipticity components e + and e x to the angular 
momentum direction L and ultimately to the tidal shear tensor by 
means of the conditional probability density p{L\<b a p)AL provides 
means to derive ellipticity correlations in terms of tidal shear cor- 
relations, and ultimately as functions of the CDM-spectrum. Be- 
cause we are interested in ellipticity correlations of galaxies, we 
smooth the CDM-spectrum with a Gaussian filter on the mass-scale 
3 x 10 u M Q //i. In this way, one can write down a 3-dimensional 
correlation function of the ellipticity field as a function of moments 
(„{r) (see Crittenden et al. 2001) of the CDM-spectrum and finally 
to carry out a Limber projection for obtaining the angular corre- 
lation function of the ellipticity field. The correlation function can 
be Fourier-transformed to yield the ii-mode and B-mode elliptic- 
ity spectra (c.f. Sect. 2.7). The related model proposed by Mackey 
et al. (2002) yields very similar shapes for the ellipticity spectra but 
predicts smaller amplitudes. 



2.5 Weak gravitational lensing 

The lensing potential if/ is defined as the line of sight projected grav- 
itational potential <I> 



Jo X 



(17) 



with the lensing-efficiency weighted galaxy distribution G(x) as a 
weighting function, 



G(x) 



(18) 



The lensing deflection angle a is obtained from the lensing poten- 
tial if/ by differentiation, a, = dflfr. Further differentiation yields the 
projected tidal field if) = djdjif/ which can be decomposed with the 
Pauli-matrices a a , because they constitute a basis for the vector 
space of 2 x 2 matrices, 



(1 -k)ct - y+o-i - y x o- 3 , 



(19) 



(c.f. Abramowitz & Stegun 1972; Arfken & Weber 2005), where 
(To denotes the 2-dimensional unit matrix. For a = 1, 2, 3, the Pauli- 
matrices have the properties cr; = <x and tr(cr a ) = 0. Due to the 
property a a ap = a- Q S a p + ie a p y a y of the Pauli-matrices, the coeffi- 
cients a a can be recovered by using a a = 4tr(i/>cr„). In particular, 
one identifies the weak lensing convergence k = |tr(i/xxo) with the 
unit matrix <x and the two components of shear y + = ^tritfxri) 
and y x = ^lr(tf!cr 3 ). The standard expression for k can be recov- 
ered with k = ^\x(ipo-a) = j Yii djdjif/ = |A fl i/r = ^divea with the 
deflection angle a = V g if/. 

The two components of lensing shear are combined to form 
the complex shear 7 = 7+ + i7x with the transformation property 
7 — > 7 exp(2i^>) under a rotation of the coordinate frame by an an- 
gle tp. Violations of the symmetry of \p are very small and might 
e.g. be caused by geodesic effects such as lens-lens-coupling and 
Bom-corrections (see Shapiro & Cooray 2006; Cooray & Hu 2002; 
Krause & Hirata 2009; Bernardeau et al. 2010, for a detailed com- 
putation). Therefore, an expansion coefficient for the contribution 
due to cr 2 which parametrises image rotation was neglected. In the 
limit of weak lensing, the galaxy ellipticity is transformed accord- 
ing to e — > 6 + 7. 

2.6 Spectrum of the weak lensing potential 

The spectrum of the weak lensing potential follows from substitut- 
ing the line of sight-expression if/ = j dxG(x)/x® mio tne Limber- 
equation (Limber 1954), 



%G(x) 2 P*,(k= e/ X ,a(x)), 



(20) 



where we used the most basic flat-sky description. The power spec- 
trum P(s,(k, a) of the gravitational potential O at the scale-factor a 
follows from the comoving Poisson equation AO = 3HiQ m /(2d)6 
and is related to the density power spectrum P s (k, a) by 



/>*(£, a) = 



3Q m D + (a)\ PsW) 



(21) 



with the Hubble distance xh = c IHq making the Ar 4 -factor di- 
mensionless. By differentiation one obtains the spectrum C a (i) = 
£ C$(f) of the deflection angle a = V e i^ and the spectrum C K (t) = 
£ 4 /4 Q,(f) of the weak lensing convergence k = ^divaa = lAgif/, 
which is equal to the ii-mode correlation function C E {£) of the 
weak lensing shear 7 in the absence of C y B {l) (c.f. the following 
section). 

The angular spectrum C^,(£) of the weak lensing potential if/ 
resulting from the Limber-projection of P®(k) is depicted in Fig. 1 
along with the spectrum C a (C) of the lensing deflection angle a and 
the weak lensing shear spectrum C^iC), for the high redshift galaxy 
sample with Zmed = 0.9 and the low redshift galaxy sample with 
Zmed = 0.3. Furthermore, predictions for the spectra using a linear 
and a nonlinear CDM-spectrum P{k) are contrasted. 



2.7 Correlations of spin-2 fields 

Both the galaxy ellipticities and the Stokes-parameters of the CMB- 
polarisation form a tensorial spin-2 field, which means that rota- 
tions of the coordinate frame by an angle ip give rise to a trans- 
formation of the tensor components as e — > exp(2ii^)e and P — > 
exp(2iif)P, when the ellipticity is written as a complex ellipticity 
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multipoie'order I ' 



Figure 1. Angular spectrum Cj,{€) (green line) of the lensing potential if/, 
the spectrum C„(t) = l'Cip({) (blue line) of the lensing deflection field 
a = Vfli/r and the spectrum C 7 E (€) = CM = Cy(f) of the £-mode 
shear, for the high redshift galaxy sample (solid lines) and the low redshift 
galaxy sample (dashed lines), comparing the prediction for a linear CDM 
spectrum (thin lines) with a nonlinear one (thick lines). 



e = e+ + ie x and the polarisation tensor P is composed of the Stokes 
parameters Q and U according to P = U + iQ. 

Correlations between two points t and 2 separated by the 
distance 6 of a spin-2 field such as the complex polarisation P = 
Q + ill, the complex ellipticity e = e+ + ie x or the weak lensing 
shear 7 = 7+ + i7 x are described in terms of two functions £ ± (9), 



f + (0) = <e + (0i)e + (0 2 )> + <e x (0 1 )e x (0 2 )>, 
= <e + (0i)f + (02)>-<e x (0i)e x (0 2 )). 



(22) 
(23) 



The correlation functions £±(6) are constructed from the variances 
of the components e + and e x using vanishing cross-correlations, 
(e + 6 x ) = 0. They can be transformed to the spectra C E {€) and C B {€) 
of the gradient (E) and vorticity (B) modes of the ellipticity field, 



cm 



cm 



n J 8d6 [& (0y o (£0) - #-(0)/ 4 (f 0)] . 



(24) 



(25) 



by Fourier transform (Kaiser 1992; Schneider et al. 2002; Schnei- 
der & Kilbinger 2007; Fu & Kilbinger 2010). Completely analo- 
gous formulae apply for the description of the angular correlation 
properties of the weak lensing shear and their transformation to 
Fourier space yielding C y E (l) and C B {€), the latter of which is zero 
if lensing on a scalar gravitational potential is considered and if the 
Born-approximation applies. 

Figs. 4 and 5 shows intrinsic ellipticity spectra C E {€) and 
C € B (£) for the EUCLID galaxy sample with its median redshift 
at Zmed = 0-9 which are contrasted with the spectra C e E (C) and 
C B (C) for a galaxy sample with a much lower median redshift of 
Zmed = 0.3. For comparison, we superpose the corresponding spec- 
tra C E {1) for the weak lensing shear 7 measured on the same galaxy 
populations, both for a linear and a nonlinear CDM-spectrum. 

The spectra are constant and equal in amplitude up to mul- 
tipoles of t - 100, indicating the absence of correlations such 
that on each scale on measures the variance of the uncorrelated 
ellipticity field. Correlations become important on angular scales 
I i 300 where the spectra level off and decrease from multipoles 



of t i 3000 on very rapidly. In the peak region, the ellipticity E- 
modes have an amplitude larger than the 6-modes by about an order 
of magnitude. Consistent with expectations the lensing spectrum 
exceeds intrinsic ellipticity correlations significantly for the high 
redshift case, but one observes the contrary behaviour for the low 
redshift case. Furthermore, the impact of nonlinear structure for- 
mation is stronger in the low-redshift case, as non-linearities had 
more time to develop. 



3 IMPRINTS ON LENSING ON ELLIPTICITIES 
3.1 Lensing effects on ellipticity fields 

In weak lensing studies it is commonly assumed that the observed 
lensed galaxies would have had no shape-correlations without lens- 
ing, and that there is no clustering in the galaxy sample, neither 
along the line-of-sight nor perpendicular to it. In short, intrinsic 
ellipticities are drawn independently from a distribution which is 
commonly assumed to be Gaussian with variance a\. In the esti- 
mation process of weak lensing spectra from galaxy shapes the in- 
trinsic shape variations would contribute the Poissonian term cr\jn, 
if the estimation process comprises n galaxies. 

The physical picture we have in mind is a background field of 
intrinsically correlated ellipticities on which gravitational lensing 
acts by deflection and shear, while there are no correlations between 
the ellipticities and the matter distribution responsible for lensing. 
In this sense, we consider lensing on II-alignments, while neglect- 
ing Gl-alignments, which vanish in the case of Gaussian fluctua- 
tions statistics for quadratic alignment models as the one used in 
this work, but would be present in nonlinear structures (Lee & Pen 
2008). Other lensing-induced effects are modulations in the surface 
density of galaxies, due to the interplay between magnification of 
the image brightness and dilution over a larger solid angle, which 
we neglect for the purpose of this paper. 



3.2 Adaptation of the CMB-lensing formalism 

As in the case of CMB lensing where the (complex) polarisation 
tensor P(6) = Q{0) + it/(0) is measured at a new position + a 
due to gravitational lensing, we assert that the ellipticity e is not 
observed at the true position of the galaxy, but at the apparent 
position 6 + a, e(ff) — » e(0 + a) with the lensing deflection angle a. 
Additionally, a variation of the deflection angle across the galaxy 
image leads to a distortion described by the complex shear 7 and 
the convergence k. 

Adapting the CMB-lensing formalism, correlations between 
the components of the shifting angle a at two positions 0i and 02 
are described by (Seljak 1996) 

<a,(0 1 )a,(0 2 )> = ic o (0) - C 2 (6) 8^ (26) 

with = 02 - 0i. The two correlation functions of the deflection 
angle are defined as 



C o (0) 



and 



r e 3 d£ 

J -^-CMM 



—cmum. 

We abbreviate the variance of the deflection angle 
0^(0) = C (0) - C o (0) 



(27) 



(28) 



(29) 
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Figure 2. Con-elation functions cr 1 (6) = Co (0) - Co (8) (blue line) and C 2 (6) 
(green line) as a function of separation angle 9, for the high redshift galaxy 
sample (solid line) and the low redshift galaxy sample (dashed line), again 
comparing the predictions from linear (thin lines) and nonlinear (thick lines) 
CDM spectra. 



in complete analogy to CMB-lensing for describing uncorrelated 
deflections. 

The characteristic function of a, i.e. the Fourier transform of 
the probability density p(a)da is then obtained as: 



(30) 



( exp (it [a(0i) - a(0 2 )]) } = exp |y [-cr 2 (8) + cos 2ip ( C 2 (6)\ 

and can be expressed in the case of Gaussian distributions in terms 
of cr 2 (Q) and Ci(#). In the case of CMB-lensing, non-Gaussian con- 
tributions have been shown to have negligible effect on the deflec- 
tion angle statistics (Carbone et al. 2009; Merkel & Schafer 2011) 
and in the case of weak cosmic shear, arguments about the rarity of 
strong deflections not described by a Gaussian distributions apply 
in a similar way (Hamana et al. 2005). 

Fig. 2 shows the quantities cr 2 (9) = Co(0) - Co(9) and C 2 i9) 
needed in this formalism, for both a high and a low redshift galaxy 
sample and for linear and nonlinear CDM-spectra: Both correlation 
functions assume larger values for the high-redshift sample and for 
lensing on nonlinear structures, which in particular causes a larger 
variance of the deflection angle on small angular scales. For large 
values of of the argument 9, both correlation function start to oscil- 
late rapidly. 

The correlation properties of the lensing-distorted ellipticity 
field can be described using the two correlation functions £±(ff), 



C(0) = (<E*(x + a)<E(x' +«')> 

g_(ff) = (exp(-4i0 f )e(x + a)e(x' + or'))- 



(3D 
(32) 



Substituting the correlation function for the deflection angle in the 
Fourier-transforms of the above expressions yields the correlation 
functions £+($) of the new ellipticity field. They can be transformed 
to £-mode and B-mode spectra with the standard transformations 
eqns. (24) and (25). 

These steps lead to a transformation formula the ii-mode and 
B-mode spectra of the ellipticity field, which can be summarised by 
a concise matrix notation: 



C' B (t) 



I™ 



W-(t, o 



W- ({,{') 



c%(n 

C B (t') 



This notation shows explicitly the mixing between scales due to the 
convolution weighted with W+(t, (') and the conversion between 
C|(f) and C%(t) under the influence of W.(t,t'). And clearly, the 
displacement mechanism can not generate ellipticity correlations 
as C' x {{) remains zero if the C x (t) is zero to begin with, X e {E, B). 
These kernels W ± (t, £') are given by 



W-(t, (') 
with the functions 

A(£, 9) = 

B{t, 6) = 



- J 9d9 [J (tff)A(t',ff) + J 4 (t0)B(t', 9)] , 

- J 0d6 [J o (t0)A(t',e) - J 4 (t0)B(t', 9)] , 



exp - 



lV 2 (6>) 
2 

t 2 o- 2 (9) 



Ut,e)+-c 2 (e)Uto) 



Mt,e)+ -c 2 ie)j s ite) 



(34) 
(35) 

(36) 
(37) 



Here, uncorrelated deflections contained in the variance cr 2 (9) give 
rise to a Gaussian convolution kernel while correlated deflections 
due to Cii9) show a more complicated mode-coupling. We abbre- 
viated J s (x) = Ji(x) + Jeix). 

In the limit absent lensing, Cq(9) = C2(#) = such that 
W+(t, f ) = S(t - t')/t and W-(t, f) = 0, due to the orthogonal- 
ity relations of the cylindrical Bessel functions, 



ite) 69 j n (te)j n (t'e) = s D (t - f )• 



(38) 



(33) 



In this case, the convolution is reduced to a Dirac So -function and 
the mixing matrix is the unit matrix, so that the £-mode and B- 
mode amplitudes are conserved and there is no convolution be- 
tween f-modes. We have verified that higher-order corrections aris- 
ing in the transformation of correlation functions do have a negligi- 
ble effect for the evolved ellipticity correlations (Challinor & Lewis 
2005; Lewis & Challinor 2006) and in our numerical implementa- 
tion, we used the same relations for the required number of grid 
points (4 x t max tabulated values of cr 2 (9) and Cii9) in 9) as sug- 
gested for CMB-lensing. 



3.3 Conversion between E and B-modes 

Fig. 3 shows the mode coupling kernels W+(t, I') and W-(t, £') for 
the high-redshift distribution and computed for a nonlinear CDM 
spectrum, with qualitatively very similar results for the low red- 
shift sample. Apart from a smooth variation of W+(t, £'), which 
acts on the ellipticity spectra by convolution, one notices tall spikes 
at £ = £', illustrating the closeness to diagonality of the W + -matrix. 
In contrast, W-(t, £') shows smaller amplitudes by about two or- 
ders of magnitude, indicating that the conversion between E- and 
B-modes is a minor effect compared to the convolution mediated by 
W+it, £'), with strong oscillatory features close to diagonal £ = £'. 
The kernels W ± (£, £') show an inverse scaling with multipole £ such 
that they become approximately constant when substituted into the 
relation (33) by multiplication with the £'d£' -differential. 



3.4 Lensed ellipticity spectra 

The final results are given in Figs. 4 and 5, which compare the 
initial ellipticity spectra C e E (£) and C e B i£) of the ellipticity field as 
predicted by correlated angular momenta, and the distorted spectra 
C' E (£) and C' B (£) that encapsulate the imprint of lensing deflection 
and therefore display altered correlation properties. For comparison 
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Figure 3. Mode coupling functions W+((,t) (green lines) and W-(t, t') 
(blue lines) used in the transformation of the ellipticity spectra, which de- 
scribes the convolving effect on the ellipticity spectra, for (' = 10, 100, 1000 
(solid, dashed and dash-dotted, respectively). The coupling functions shown 
are the ones for the high redshift galaxy sample, while those for the low red- 
shift sample look qualitatively very similar. 



Figure 4. Ellipticity spectra C%(t) (blue line) and C B (C) (green line) as 
predicted by the angular momentum model (dashed lines), and the Iensed 
ellipticity spectra (solid lines) where the deflections were computed for the 
high redshift galaxies. For comparison, we plot the spectrum C E (C) of the 
weak lensing shear y (red line). 



with weak lensing, we plot the weak convergence spectrum C K (C) 
expected from the EUCLID galaxy sample in comparison, for a 
nonlinear CDM spectrum (using the parametrisation by Smith et al. 
2003). 

The first observation is that ellipticity correlations reach am- 
plitudes similar to those of the weak lensing convergence in the 
nonlinear part corresponding to amplitudes I i 300, and that the 
intrinsic £-mode spectrum C E {€) is larger than the S-mode spec- 
trum Cg(f) by about an order of magnitude in this regime. On larger 
angular scales, there are no appreciable ellipticity correlations and 
one effectively observes the variance of the ellipticity field for un- 
correlated objects. Consequently, the spectra have identical ampli- 
tudes and are effectively constant. In this regime, the shifting ef- 
fect is not able to affect the galaxies, which is a well-known re- 
sult in CMB-lensing, where scale free-spectra are invariant (Lewis 
& Challinor 2006): The mode-conversion mechanism is ineffective 
if the spectra are equal, C E {€) = C%(£), and the convolution with 
W+({, {') is not able to redistribute amplitudes. In contrast, both 
spectra are affected on multipoles I > 1000, where in particular 
C' B (C) has decreased relative to C„{€). 

Fig. 6 shows the changes in the spectra as a function of mul- 
tipole C by giving the ratio of the evolved and initial ii-mode and 
S-mode spectra, C E {i)IC\(l) and C' B {€)IC%{€) respectively. As al- 
ready indicated by Figs. 4 and 5, we see a significant decrease of 
amplitude amounting to 5% for the E- and 30% for the B-modes 
from t = 3000 on in the case of the high-redshift sample and from 
( = 1000 on in the case of the low-redshift sample. This implies that 
for EUCLID's weak lensing application, changes in the ellipticity 
spectra are affecting scales where the shape noise starts dominat- 
ing, but for shallower surveys, lower multipoles would be affected 
by weak lensing deflection. We conclude that in the case of deep 
surveys such as EUCLID, weak lensing manifests itself primarily 
as weak lensing shear which dominates over intrinsic alignments 
and the lensing deflection effect shapes intrinsic alignments by de- 
creasing their amplitudes only at very high multipoles. 




C*(l), linear 

C E (I), Iensed 

G B (I), Iensed 



3 10 30 J 00 300 1000 3000 

multipole order t 

Figure 5. Ellipticity spectra C E (C) (blue line) and C B (f) (green line) as 
predicted by the angular momentum model (dashed lines), and the Iensed 
ellipticity spectra (solid lines) where the deflections were computed for the 
low redshift galaxies. For comparison, we plot the spectrum C £ (() of the 
weak lensing shear y for the same galaxy sample (red line). 



3.5 Violated homogeneity of the ellipticity field 

Lensing of the intrinsic ellipticity field introduces a violation of 
their statistical homogeneity in complete analogy to the lensing de- 
flection acting on the CMB polarisation. For a given realisation of 
the deflection potential 41 one can estimate the lensing effect on the 
ellipticity spectra, X = E,B: 

w x (tK(t')) = fx((, tmt - o (39) 

with 

f E {t,t') = (t-e')[tc e E (€) + e'c E (e)]cos2<p e , r m 

f B (e J') = {t-t')[tC c B (t)+t'C\(t')\cos2y (fi (41) 
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Figure 6. Ratios C' E (C) I C € E {C) (blue lines) and C' B (e)/C B ({) (green lines), 
for both a high redshift galaxy sample (solid lines) and a low redshift galaxy 
sample (dashed lines). 

where tpc.c the enclosed angle. Nonzero correlations between mul- 
tipoles are the signature of homogeneity violation introduced by a 
single realisation and they would disappear in the process of en- 
semble averaging the lensing potential. 

In order to place an upper limit on this effect we select a par- 
ticularly simple geometry, namely parallel alignment of the wave 
vectors t and {' such that the cosines are equal to one and have the 
wave vectors them differ by one unit, t — (' = 1, as the coupling 
between neighbouring multipoles is strongest due to the rapid de- 
cline of the lensing potential's Fourier transform with increasing I. 
In this limit, and if one assumes the ellipticity spectra to be slowly 
varying, the ratio between the off-diagonal and diagonal correla- 
tions is given by 

- E — = - J — = —;crAl = 1), (42) 

C\(f) C\{€) 2tt 2 * 

where we replaced tp by cr^, as an order of magnitude estimate for 
a typical amplitude. Substituting numbers yields upper limits of 
10% for the off-diagonal correlation relative to the diagonal ones 
at t = 10 3 which is smaller than the weak lensing shear at high 
redshifts, but comparable if not slightly larger than the lensing shear 
at low redshifts. 

As the imprints of weak lensing on the intrinsic ellipticity pat- 
tern are basically identical to those in the case of lensing of the 
CMB polarisation, it is conceivable that the lensing deflection field 
can be estimated by measuring the amount of off-diagonal (mean- 
ing I -t (') spectra of the ellipticity modes as illustrated above, 
by applying the reconstruction technique worked out by Hu and 
Okamoto: With a model for intrinsic ellipticity spectra (which can 
be well predicted using a good prior on Q.,„, <x 8 , a and a) the statis- 
tics of the lensing deflection field can be inferred from broken ho- 
mogeneity. In contrast to polarisation tensors, however, the ellip- 
ticity field is strongly shaped by shear and only in extreme cases 
such as a low-redshift galaxy sample at large multipoles the intrin- 
sic alignment effects are dominant, where of course issues with the 
low surface density of lensing galaxies and the corresponding high 
Poisson-noise become important. 

For illustration, Fig. 7 shows a realisation of the lensing po- 
tential with derived lensing deflections as gradients of the poten- 
tial. The second derivatives have been used for generating a shear 




-6 -4 -2 2 4 6 



Figure 7. Realisation of a lensing potential on a 14° X 14° -patch with 100 
randomly sampled galaxies. The ellipticity is given by the lensing shear and 
the arrow indicates the apparent shift due to the weak lensing deflection. 
The realisation uses a lensing potential for the high-redshift galaxy sample 
with dominating shear. The displacement arrows are enhanced by a factor 
of 100 and the shears by a factor of 10, applied to a circular spot. 

field that is depicted as a shape distortions of the otherwise circu- 
lar spots. For visualisation, the deflections have been enlarged by a 
factor of 100 and the shear has been multiplied by 10. 



4 SUMMARY 

The subject of this work was an investigation of lensing effects and 
their observable signatures if the lensed galaxy sample shows in- 
trinsic ellipticity correlations. Apart from weak lensing shear that 
operates on the shape of galaxies there will be a lensing deflection, 
which is unobservable in the case of uncorrelated ellipticities, but 
generates observable signatures if the lensing galaxies are intrin- 
sically shape correlated. The lensing deflection manifests itself in 
the ellipticity spectra in three distinct ways: Firstly, there is a loss of 
amplitude in the spectra at high multipoles, secondly one observes a 
redistribution of amplitude between the E-mode and S-mode spec- 
tra, and thirdly there will be correlations between adjacent multi- 
poles. In deriving these effects we made heavy use of analogies 
to the theory of lensing of the cosmic microwave background po- 
larisation and identical mathematical properties of the ellipticity- 
the polarisation tensors. In our investigation, we are comparing the 
forecasts for a high-redshift lensing survey such as EUCLID with 
a low-redshift galaxy sample in order understand the scaling be- 
haviour of all effects with distance. 

(i) We derive ellipticity ii-mode and B-mode spectra from a 
physical alignment model due to Crittenden et al. (2001). This 
model is quadratic in the tidal shear and applicable for describing 
shape correlations between isolated spiral galaxies. Due to the lack 
of an analytical description of how a galactic disk is oriented inside 
a dark matter halo we think of our spectra as upper limits as we 
assume perfect alignment of the symmetry axis of the galactic disk 
with the host halo's angular momentum direction. The two parame- 
ters that enter our ellipticity model, the alignment parameter a and 
the disk thickness a, are determined from numerical simulations 
and from observations of local galaxies. 
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(ii) The impact of lensing deflection on intrinsically shape- 
correlated galaxies is threefold: There is a smoothing of the in- 
trinsic ellipticity spectra, a mixing in multipole and a conversion 
between E-modes and S-modes, and the generation of correlations 
between otherwise uncorrelated multipoles, as an expression of vi- 
olated homogeneity of the lensed galaxy field. 

(iii) By drawing analogies between galaxy ellipticities and the 
CMB-polarisation, namely that both are tensorial fields with spin-2, 
we can formulate transformation formulas for the ellipticity spec- 
tra, if individual galaxies have been coherently shifted to a new po- 
sition by lensing. The transformation formula can be written con- 
cisely as a combined convolution and mode-mixing relation. 

(iv) Lensing deflection operates on intrinsic ellipticity spectra 
by convolution. Correlation amplitudes are redistributed in multi- 
poles which can be observed on small angular scales when intrin- 
sic alignment spectra cease to be constant and drop in amplitude. 
Then, lensing causes the spectra to drop faster. Qualitatively, these 
effects are weak at high redshifts and dominated by far by the weak 
shear signal, but are sizable at low redshifts, where the weak lens- 
ing shear is small. In this case, weak lensing can actually weaken 
shape-correlations by random redistribution of intrinsically aligned 
galaxies. 

(v) The losses in amplitude amount roughly to 5% in C e E ({) and 
to 30% in C £ B (£) at £ = 3000 in the case of the high redshift sample 
and at I = 1000 in the case of the low redshift sample. Compared 
to the convolution of the spectra the conversion between E- and 
S-modes is a minor effect. 

(vi) We have derived an upper limit on the correlation between 
different multipoles due to broken homogeneity by using the fact 
that the correlations with between adjacent multipoles should be 
strongest. These correlations can be estimated to be at most ~ 10% 
of the spectra at the largest multipoles, both for E- and S-modes, 
and are proportional to I. 

(vii) Although we could take advantage of formal analogies be- 
tween the CMB-polarisation and ellipticity fields, concerning sym- 
metry properties, the description with spectra and the incoiporation 
of the lensing effect we would like to emphasise that that in contrast 
to the CMB, lensing does not introduce a bispectrum into the ellip- 
ticity correlations. The CMB-lensing bispectrum is sourced by the 
integrated Sachs-Wolfe effect in the same potential that causes the 
lensing-deflection (Hu 2000), and there is no analogous mechanism 
in the case of galaxy ellipticities. 

We conclude that for deep-reaching lensing surveys intrinsic align- 
ments are subdominant and that the shaping of their correlations 
by weak lensing deflection (and by peculiar motion, which is of 
a similar order of magnitude, see Giahi-Saravani & Schafer 2012) 
is small compared to gravitational shear. At low redshifts, how- 
ever, the situation is inverted: Intrinsic alignments dominate and 
the most important lensing effect is deflection. In contrast to the 
CMB-polarisation, it is doubtful if a violation of homogeneity of 
the ellipticity field introduced by lensing can be observed. 
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